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1. Introduction 

The observations of neutrino oscillation from the sun, the atmosphere, reactors, and particle 
beams have established that neutrinos have mass and that three flavor states (v^, v^-) are super¬ 

positions of three mass states (Vi, V 2 , V 2 ). The mixing can be quantified using a 3 x 3 matrix, known 
as the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix, which can be parameterized in 
terms of three mixing angles {dn, 623 , ^ 13 ) and a CP violating phase (5). The neutrino oscillation 
also depends on two mass-squared differences Arn^i and Amjj. 

Among the three mixing angles, the value of 613 can be extracted using the electron antineu¬ 
trinos from the reactors at a short baseline (~km). The survival probability is given by 

Pve-^Ve = 1 — cos^6i3sin^26i2sin^A2i — sin^20i3(cos^6i2sin^A3i-|-sin^6i2sin^A32), (1.1) 

where Ay/ = 1.261 and Amji is the difference between the mass-squares of the 

mass eigenstates Vy and v,. Besides the mixing angle 613 , the effective mass-squared difference, 
Aml^, defined by sin^ Age = cos^ 612 sin^ A 31 -|- sin^ dn sin^ A 32 , can also be determined by measure¬ 
ment of the energy-dependent oscillation. 

In 2012, the Daya Bay experiment observed a non zero 613 at a significance of 5.2a [1] via a 
rate-only measurement with 55 days of data. Other experiments observed consistent results [2, 3,4, 
5]. In 2013, Daya Bay updated the results with 217 days of data, with a spectral analysis to measure 
the oscillation frequency, which led to the first direct measurement of the effective mass-squared 
difference Am^g in the electron antineutrino disappearance channel. The relative large value of 
013 and its precise measurement facilitate future neutrino oscillation research on the neutrino mass 
hierarchy and CP violating phase. 

2. The Daya Bay Experiment 

The Daya Bay Reactor Neutrino Experiment was designed to make a precise measurement 
of the mixing angle ^ 13 . The precision is ensured by the following features: 1) high statistical 
precision due to large thermal power of reactors and large target mass of detectors; 2 ) reduction of 
reactor related uncertainties by relative Far/Near measurements; 3) reduction of detector related un¬ 
certainties by identically designed multiple detectors; 4) low background by sufficient overburden 
and good shielding. 

A detailed description of Daya Bay experiment can be found in Refs. [ 6 , 7]. The Daya Bay 
experiment is located near three nuclear power plants (Daya Bay, Ling Ao, and Ling Ao II) on the 
southern coast of China. Each power plant has a pair of reactor cores with a distance of 88 m. All 
six cores are functionally identical pressurized water reactors, each with a maximum of 2.9 GW 
thermal power. A total of 8 identically designed antineutrino detectors (ADs), with a 20-ton target 
mass each, are placed in three underground experimental halls (EHl and EH2 for the near sites, 
and EH3 for the far site), covering baselines ranging from 360 m (EHl to Daya Bay) to 1910 m 
(EH3 to Daya Bay). 

Each AD consists of three nested cylindrical zones separated by two concentric acrylic vessels. 
The innermost zone contains 20 tons of gadolinium-doped liquid scintillator (Gd-LS). The middle 
zone is filled with 21 tons of un-doped liquid scintillator (LS), while the outer zone consists of 37 
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tons of mineral oil. The Vg is detected via the inverse j3-decay (IBD) reaction, Vg + /r —)• + n, 
in the Gd-LS. The coincidence of the prompt positron light with the delayed gamma rays (totalling 
~8 MeV) generated from the neutron capture on Gd with a mean capture time of ~30 /is provides 
a distinctive Vg signature. In each AD, light created as a result of particle interactions in the Gd- 
LS and LS, is collected by 192 radially positioned 20-cm PMTs in the MO. The detectors have a 
light yield of ~165 photoelectrons/MeV and a reconstructed energy resolution of ss 8% at 
1 MeV. Three automated calibration units (ACUs) mounted on the stainless steel vessel lid allow 
for remote deployment of a light-emitting diode (LEDs) and calibration sources, such as ^^Ge, 
^‘^^Am^^C and ^°Co, into the Gd-LS and LS volumes along three vertical axes. The ADs in each 
EH are shielded with > 2.5 m of high-purity water against ambient radiation in all directions. 
Each water pool is segmented into inner and outer water shields (IWS and OWS) and instrumented 
with photomultiplier tubes (PMTs) to function as Cherenkov-radiation detectors whose data were 
used by offline software to remove spallation neutrons and other cosmogenic backgrounds. The 
detection efficiency for long-track muons is >99.7% [8]. The water pool is covered with an array 
of resistive plate chambers (RPCs). 

The detector energy scale of all ADs was calibrated using ^^^Am*^C sources providing an 
~8 MeV peaks from neutron capture on Gd. The time variation and the position dependence of 
the energy scale was corrected using the 2.506 MeV gamma-ray peak from ^‘’Co sources. The 
reconstructed energies of various calibration sources in different ADs are compared in Pig. 1. The 


o Neutron from muon spallation a Alpha from natural radioactivity 

• Neutron from IBD ■ Gamma from calibration source 

O Neutron from Am-C source ♦ Gamma from natural radioactivity 



Figure 1: Comparison of the reconstracted energy of calibration sources between antineutrino detectors. 
Bad is the reconstructed energy determined for each AD, and (E) is the 8-detector average. Error bars are 
statistical only, and the difference between detectors for all calibration sources are < 0.2%. 


ACU calibration sources Am^^C and ^‘^Co were positioned at the detector center. Neutrons from 
IBD and muon spallation captured on Gd were nearly uniformly distributed in the Gd-LS region. 
a particles from polonium decays with positions estimated to be within the Gd-LS region provided 
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another diffuse reference. Comparisons of intrinsic and spallation neutron capture on 

hydrogen were biased by backgrounds, therefore only those interactions with reconstructed posi¬ 
tions within 1 m of the center of each detector were considered. Based on all this information, 
the uncorrelated uncertainty of the energy scale was determined to be 0.2% in the energy range of 
reactor Vg. 

The energy response of the ADs is not linear due to scintillator and electronics effects, each 
contributing at the level of 10%. The scintillator nonlinearity is particle and energy dependent, and 
is related to intrinsic scintillator quenching and Cherenkov light emission. The electronics non¬ 
linearity is introduced due to the interaction of the scintillation light time profile and the charge 
collection of the front-end electronics. The energy nonlinearity model was determined by a com¬ 
bined fit to monoenergetic y peaks from radioactive sources and the continuous y + p spectrum 
extracted from produced by a spallation neutron after muon interactions with carbon. The 
electronics nonlinearity was determined by studying the time profile of charge in fhe dafa and MC. 
The uncerfainfies of fhe nonlinearify are < 1 % above 2 MeV. 

The Daya Bay experimenf sfarfed dafa-faking on September 23, 2011 wifh fwo ADs insfalled 
in EHl. On December 24, 2011, fhe "6-AD" run sfarfed shorfly affer four more ADs were insfalled 
in EH2 and EH3. The previous resulfs [1, 9, 10] on fhe measuremenfs of oscillafion parameters are 
based on fhe 6-AD dafa. The "8-AD" run sfarfed on October, 2012 after fhe insfallafion of fhe final 
fwo ADs in fhe summer of 2012. The 8-AD dafa obfained up to Nov. 27, 2013 (corresponding 
to 404 days) were analysed fogefher wifh fhe 6-AD dafa fo exfracf fhe lafesf resulfs of 0i3 and 
A/Wg^ [11]. An independenf measuremenf of On using fhe neufron capfure on hydrogen of fhe 
6 -AD dafa has also been performed [12]. The measuremenf of fhe reactor anfineufrino flux and 
specfrum has been obfained for fhe 6-AD dafa [13]. Moreover, a search for lighf sterile neufrino 
using fhe 6-AD dafa has been reported [14]. In fhe nexf secfion, we will describe fhese recenf 
resulfs. 

3. Recent Results 

3.1 Oscillation Analysis from Neutron Capture on Gadolinium 

A fofal of ~ 1.1 million (~150 k) IBD candidates for fhe near (far) sites were selected using fhe 
same criferia described in [1]. The IBD rate as well as fhe prompf and delayed energy specfra show 
consisfency befween side-by-side detectors. The relative energy scale uncerfainly was improved fo 
0.2%. This reduction of 43% compared fo previous publication [10] was achieved by improvemenfs 
in fhe correcfion of posifion and lime dependence, and enhanced fhe precision of Aml^ by 9%. 

Eslimales for fhe five major sources of background for fhe new dafa sample are improved wifh 
respecf to Ref. [10]. Two of fhe fhr'ee Am^^C sources in each AD in EH3 were removed during 
fhe 2012 summer insfallafion period. As a resull, fhe average correlaled ^"^^Am^^C background 
rale in fhe far hall decreased by a facfor of 4 in fhe 8-AD period. The eslimafe of fhe fasl neufron 
background was improved by fagging fhe fasl neufron candidales following cosmogenic signals 
defected by fhe OWS or RPC. The energy specfrum of fhese velo-lagged signals was consislenl 
wifh fhe specfrum of IBD-like candidale signals above 12 MeV, and was used fo eslimafe fhe rate 
and energy specfrum for fhe fasl neufron background in fhe range of 0.7-12 MeV. The melhods in 
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Refs. [1, 9] were used to estimate the backgrounds from the accidental background, the correlated 
— n decays from cosmogenic ^Li and ^He, and the ^^C(a, n)*^0 reaction for the current 6+8 AD 
data. The total backgrounds amounted to about 3% (2%) of the IBD candidates in the far (near) 
hall(s). 

The reconstructed prompt energy spectrum observed in the far site is compared with the ex¬ 
pectation based on the near-site measurements in left plot of Fig. 2. In the right plot of Fig. 2, 
the measured survival probability as a function of L/E from all experimental halls is compared 
to the expectation assuming the best fit oscillation parameters. The spectrum distortion is highly 
consistent with the oscillation interpretation in the three neutrino framework. 




Figure 2: Left; (Top panel) Background-subtracted reconstructed positron energy spectrum observed in the 
far site (black points), as well as the expectation derived from the near sites excluding (blue line) or including 
(red line) the best estimate of oscillation. (Bottom panel) Ratio of the spectra to the no-oscillation case. The 
error bars show the statistical uncertainty of the far site data. The shaded area includes the systematic 
and statistical uncertainties from the near site data. Right: Electron antineutrino survival probability versus 
effective propagation distance Leff divided by the average antineutrino energy {Ef). The data points represent 
the ratios of the observed antineutrino spectra to the expectation assuming no oscillation. The solid line 
represents the expectation using the best estimates of sin^20i3 and \A.m^f\. The error bars are statistical only. 


Oscillation parameters are extracted by a relative measurement of prompt energy spectra at 
the near and far halls. A covariance matrix including known systematic and statistical uncertainties 
was used in a fit. We found sin^20i3 = 0.084±0.005 and = 2.421q}j x 10^^ eV^ with 
Z^/NDF = 134.7/146. Fig. 3 shows the allowed regions for sin^20i3 and lA/n^^l at the 68.4%, 
95.5%, and 99.7% confidence levels. If is worth noting fhaf fhe currenf accuracy is sfill dominafed 
by sfafisfics af fhe far site. The jAmg^l measuremenf is highly consisfenf wifh and of comparable 
precision fo fhe muon neufrino disappearance experimenfs [15, 16]. Under fhe normal (inverted) 
hierarchy assumption, |Amgg| yields |A;m 32 | = 2.39lQjj eV^ (|Am 32 | = —2.49+q}j eV^). 

3.2 Oscillation Analysis from Neutron Capture on Hydrogen 

An independent measurement of sin^20i3 has been obtained via the detection of IBDs tagged 
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Figure 3: Regions of 68.3%, 95.5% and 99.7% confidence levels in the |AOTgg|-sin^20i3 parameter space. 
The adjoining panels show the dependence of on sin^20i3 (top) and |Amgg| (right). The |Amgg| allowed 
region (shaded band, 68.3% C.L.) was consistent with measurements of |Am 32 | using muon neutrino disap¬ 
pearance by the MINOS [15] and T2K [16] experiments, converted to |AOTgg| assuming the normal (solid) 
and inverted (dashed) mass hierarchy. 


by neutron capture on hydrogen (nH). Comparable statistics as the nGd case is achieved by the 
~15% of neutron captures in the Gd-LS region and almost all of the neutron captures in the LS 
region. New analysis approaches have been developed to meet the challenges associated with the 
higher background, longer neutron capture time (~200 /is), and a lower energy (2.2 MeV) 7 ray 
from neutron capture for nH IBD events. The prompt candidate is required to have £ > 1.5 MeV, 
a longer coincidence time window of 400 /is is required, and the distance between the prompt 
and delayed candidates is required to be < 0.5 m. Based on the rate deficit observed at the far-site 
detectors for the 6 -AD period, a value of sin^26i3 = 0.083 ±0.018 is extracted. The result is shown 
in Fig. 4, which is consistent with the result obtained in the nGd analysis. 

3.3 Measurement of the Reactor Antineutrino Flux and Spectrum 

In the complete 6 -AD period, 296721 (41589) IBD candidates were detected in the near (fall) 
halls, and were used for a precise measurement of reactor antineutrino flux and spectrum. The 
measured IBD yield, expressed as T = (1.55 ±0.04) x 10^^^ cm^/GW/day, or Of = (5.92±0.14) x 
10^^^ cm^/fission, is shown in the left top plot of Fig. 5. The measured IBD yields are consistent 
among all ADs. A correlated uncertainty in the detection efficiency of 2.1% is the dominant factor 
for the overall uncertainty of 2.3%. A global fit for R, defined as the ratio of the Daya Bay measured 
rate to the model prediction, was performed to compare with the previous reactor antineutrino 
flux measurements following the method described in Ref. [17]. The Huber [18] and ILL [19, 
20] models predict the reactor Vg spectra for ^^^Pu, and ^^*Pu, while the Mueller [21] and 
Vogel [22] models predict it for We obtained R = 0.946 ± 0.022 (0.991 ± 0.023) for the 
Huber-i-Mueller (ILL-i-Vogel) model, shown in the left bottom plot of Fig. 5. 

In addition to the reactor antineutrino flux, the energy spectrum was measured and compared 
with the model prediction shown in the right plot of Fig. 5. The spectral discrepancy around 5 
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Figure 4: The observed prompt energy spectrum of the far hall ADs (blue) and near hall ADs (open circle) 
weighted according to baseline. The far-to-near ratio (solid dot) with best ht 613 value is shown in the lower 
panel. In the inset is the ratio of the measured to the predicted rates in each AD vi. baseline, in which the 
baselines of far ADs were shifted relative to each other by 30 (—30) m for clarity. 

MeV prompt energy is visible. The discrepancy is found to be time-independent and correlated 
with reactor power, therefore disfavoring hypotheses involving detector response and unknown 
backgrounds. A recent ab-initio calculation of the antineutrino spectrum shows a similar deviation 
in the 4-6 MeV energy region [23]. 

3.4 Search for a Light Sterile Neutrino 

A search for light sterile neutrino mixing was performed in the 6-AD period. If the sterile 
(fourth) neutrino exists, its presence could be detected via the modification to the three-neutrino 
oscillatory behavior. The experiment’s unique configuration of multiple baselines from six reactors 
to six ADs makes it possible to test for oscillations to a sterile neutrino in the 10^^ eV^ < | Anj|j | < 
0.3 eV^ range. The relative spectral distortion due to electron antineutrino disappearance was found 
to be consistent with that of the three-flavor oscillation model. Fig. 6 shows the derived limits on 
sin^ 20 i 4 covering the 10^^ eV^ < lAnr^j | < 0.1 eV^ region, which was largely unexplored. 

4. Summary 

After the final two of eight antineutrino detectors were installed in the summer of 2012, an 
additional 404 days of data was acquired with the full 8-AD configuration. The uncertainties of 
sin^ 20 i 3 and |Am^g| were halved as a result of improvements in calibration, background estimation, 
as well as increased statistics. An analysis of the relative antineutrino rates and energy spectra 
between detectors yielded sin^20i3 = 0.084±0.005 and \Amlg\ = 2.42+q}j x 10^^ eV^. This is 


























Recent Results from Daya Bay 


Liang Zhan 



10'I 
10 "S 
10 -=^- 
10^1 
lO-®-" 
10 -® 


Figure 5; Left: (Top panel) The measured Vg rate in the six ADs (closed circles), with correction of flux- 
weighted fission fractions at different sites (open squares). The average ( uncertainty) of three near detectors 
is shown in gray line (band). The rate predicted by Huber-rMueller (ILLH-Vogel) model and its uncertainty 
are shown in blue (orange) region. (Bottom panel) The measured reactor Vg rate as a function of the distance 
from the reactor, normalized to the Huber-rMueller model. The blue shaded region represents the global 
average and its 1(7 uncertainty. The 2.7% model uncertainty is shown as a band around unity. Right: (Panel 
A) Predicted and measured prompt energy spectra. The prediction is based on the Huber-rMueller model. 
The gray hatched (only in panel B) and red hlled bands represent the square-root of diagonal elements of 
the covariance matrix (\/ (Vu)) for the reactor related and the full systematic uncertainties, respectively. The 
error bars on the data points represent the statistical uncertainty. (Panel B) Ratio of the measured prompt 
energy spectrum to the predicted spectrum (Huber-rMueller model). The blue curve shows the ratio of the 
prediction of ILL-rVogel model to the Huber-rMueller model. (Panel C) The defined distribution (Xi) of 
each bin (black solid curve) and local p-values for 1 MeV energy windows (magenta dashed curve). See 
Ref. [13] for the definitions of those variables. 

the most precise measurement of sin^ 20 i 3 to date and the most precise measurement of of |Anigg| 
via electron antineutrino disappearance. The Daya Bay experiment is scheduled to continue data- 
taking with all 8 ADs until the end of 2017 at least, when the precision on both parameters is 
expected to reach 3%. 

Several other analysis results are presented. An independent measurement of sin^20i3 via 
the detection of IBDs tagged by neutron capture on hydrogen yielded sin^26i3 = 0.083 ±0.018, 
which was consistent with the nGd result. A precise measurement of reactor antineutrino flux 
and spectrum has been performed. A search for sterile neutrino mixing set stringent limits in the 
10^^ eV^ < |Am 4 j | < 0.3 eV^ range, which was largely unexplored. An improvement in the reach 
of all of these analyses is expected as more statistics are collected and the systematic errors is 
reduced. 
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Figure 6: Exclusion contours for the neutrino oscillation parameters sin^20i4 and Normal mass 

hierarchy is assumed for both Amj j and Am^j. The red long-dashed curve represents the 95% C.L. exclusion 
contour with the Feldman-Cousins method [24]. The black solid curve represents the 95% CLs exclusion 
contour [25, 26]. The parameter-space to the right side of the contours are excluded. For comparison, 
Bugey’s [27] 90% C.L. limit on Ve disappearance is also shown as the green dashed curve. 
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